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We discuss the ground state of the spin-orbital model for spin-one ions with partially filled t2g 
levels on a honeycomb lattice. We find that the orbital degrees of freedom induce a spontaneous 
dimerization of spins and drive them into nonmagnetic manifold spanned by hard-core dimer (spin- 
singlet) coverings of the lattice. The cooperative "dimer Jahn- Teller" effect is introduced through 
a magnetoelastic coupling and is shown to lift the orientational degeneracy of dimers leading to a 
peculiar valence bond crystal pattern. The present theory provides a theoretical explanation of non- 
magnetic dimerized superstructure experimentally seen in Li2Ru03 compound at low temperatures. 
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The nearest-neighbor Heisenberg antiferromagnet on a 
bipartite lattice has a Neel-type magnetically long-range 
ordered ground state. However, such a classical order of 
spins can be destabilized by introducing a frustration into 
the system through the competing interactions that may 
lead to the extensively degenerate classical ground states 
[H . In such systems exotic quantum phases without long- 
range order can emerge as the true ground states. In 
this Letter we want to point out and discuss another 
scenario, that can appear when magnetic ions on a bi- 
partite lattice possess also an orbital degeneracy. The 
physics of such systems may be drastically different from 
that of pure spin models, as the occurrence of an orbital 
ordering can modulate the spin exchange and preclude 
the formation of magnetically ordered state on a bipar- 
tite lattice. In the following we focus on a system with 
threefold-orbitally-degenerate 5=1 magnetic ions on a 
honeycomb lattice. This model is suitable to describe d? 
and -type transition-metal compounds with partially 
filled t2g levels, like the layered compound Li2Ru03 [2]. 
Here the layers are formed by edge-sharing network of 
RuOe and LiOe octahedra. The Ru ions make a honey- 
comb lattice and Li ions reside in the centers of hexagons. 
These layers are well separated by the remaining Li ions. 
The magnetically active Ru^+-ions are characterized by 
four electrons in the threefold degenerate i2g-manifold 
coupled into a S" = 1 state. Li2Ru03 undergoes a metal- 
to-insulator transition on cooling below 540 K Q . At the 
transition the magnetic susceptibility shows a steep de- 
crease and its low temperature value can be considered 
to be due almost entirely to the Van Vleck paramag- 
netism. The structural analyses have revealed the forma- 
tion of dimerized superstructure of Ru-Ru bonds in the 
low temperature phase. These observations indicate that 
Ruthenium spin-one degrees of freedom are mysteriously 
missing at low temperatures and suggest the formation of 
an unusual spin-singlet dimer phase in the ground state 
of the system. 

Here we describe the microscopic theory behind the 



stabilization of such a spin-singlet dimer state. We argue 
that, remarkably, such a novel phase can be realized on a 
honeycomb lattice because of orbital degeneracy, without 
invoking any exotic spin-only interactions. A possibility 
of formation of orbitally driven magnetically disordered 
states has been strgg^ted within various coupled spin- 
orbital models 0, Ull, @, S H, i]- The orbital induced 
frustration in t^g based systems on a bipartite (cubic) lat- 
tice has also been considered [13, The emergence of 
new phases due to the p-orbital degeneracy of cold atoms 
in optical lattices has been recently discussed within a 
spinless fermion model on a honeycomb and other two- 
dimensional lattices [13, E!]. Yet, the peculiar case of 
partially filled t2g levels on a honeycomb lattice leads to 
new results: the onset of an orbitally driven spin-singlet 
dimer phase in a spin-one system. 

The model - We assume that the low-temperature insu- 
lating phase of Li2Ru03 is of Mott-Hubbard type and de- 
scribe the low energy physics within the Kugel-Khomskii 



type spin-orbital Hamiltonian 1J|. We consider undis- 
torted honeycomb lattice of Ru ions and look for pos- 
sible instabilities towards symmetry reductions. In the 
Li2Ru03 crystal structure three distinct bonds of hon- 
eycomb lattice are in xy, xz, and yz planes (in cubic 
notations). We consider the leading part of the nearest- 
neighbor (NN) hopping integral of t2g orbitals {d^y , dy^ , 
and dxz) due to the direct cr-type overlap. The dda over- 
lap in aP plane (a, (3 = x, y, z) connects only the or- 
bitals of same ajd type. The effective spin-orbital Hamil 
Ionian for such a system been reported in Refs. 
It has the following form: 



H = 



^ • S, - l]Oy - [JoS, ■ S, + Ji]Oy }, (1) 



where the sum is taken over pairs of NN sites, Si are 
spin-one operators, and the orbital contribution are de- 
scribed by Oij and Oij operators. The second-order vir- 
tual processes locally conserve orbital index. The orbital 
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degrees are thus static Potts-like variables and their con- 
tribution can be expressed simply in terms of projectors 
Pi.Q^ onto the singly occupied orbital state a/? at site i. 
With this definition of the projectors the orbital part of 
the Hamiltonian can be written in the form equally valid 
for and its particle-hole symmetry related config- 
urations. The orbital operators Oy and Oy along the 
bond ij in a/3-plane are given by: Oij — Pi^af3Pj,ai3 and 
Oij = Pi,Q/3(l - P'jMf}) + P]Mi3{l - Pi,ai3)- For further 
analysis it is convenient to rewrite the Hamiltonian as 
the sum of three terms: H ~ Eq + -Haf + -^^fm, 

-H^AF =y~^ J[Si ■ Sj + C] Oij , HpM =—Jo ^ Si ■ SjOij , 

(ij) (u) 

(2) 

Eq = — 2JiiV, N is number of lattice sites, and C, — 
2Ji/J — 1 [17]. For the full expressions of the exchange 
constants in terms of hopping integral i, on-site Coulomb 
repulsion U and Hund's coupling Jh we refer reader 
to Refs. [H, To the leading order in small pa- 

rameter rj ~ Jb/U, the coupling constants are given 
by J » (1 - 7y)^, Jo « r?|^, Ji « (1 + 2r;)^, and 
C ~ 1 -I- 677. For our further analysis, it will be important 
that C > 1 for all values of the Hund's coupling. The 
inequality implies that any magnetically ordered ground 
state has a positive classical energy on antiferromagnetic 
(AF) bonds. It suggests that the formation of a low di- 
mensional network of AF spin-coupling patterns are ener- 
getically favorable due to a larger gain of quantum spin- 
energy per bond. This picture is substantiated by the ex- 
act solution of the ground state problem at zero Hund's 
coupling described in the next section. 

Zero Hund's coupling - The spin-orbital model defined 
in Eq. ([2]) has a small parameter <C 1. We start our 
analysis from the limit of zero Hund's coupling 7/ = 
(_ffFM=0) and look for the possible spin-coupling patterns 
generated by Haf term of the Hamiltonian. 

The antiferromagnetic term Haf is active only on the 
bonds with corresponding orbital being singly occupied 
at both ends of a bond (e.g on a bond ij in xy plane 
both sites have xy orbital singly occupied). The AF 
bonds can only form non-intersecting linear open and/or 
closed chains (see Fig. [TJ\). The longer is the chain, the 
more energy can be gained from the AF spin interaction. 
However, for each AF bond we pay a positive energy 
( > I (the second term in Haf)- Under these condi- 
tions the minimal possible AF energy is achieved when 
all AF chains are dimers. The proof of the stability of 
dimer states against the formation of longer open Heisen- 
berg chains is based on the variational estimate on the 
ground-state energy of the M-site spin-one AF Heisen- 
berg chain with open ends: Em > 1 — with the 
equality attained only at M = 2. This estimate can be 
obtained by dividing the chain into shorter overlapping 
sub-chains of lengths two and three with exactly known 




FIG. 1: (Color online) Examples of orbital and spin-coupling 
patterns on the honeycomb lattice of Ru ions. A) Decoupled 
AF chain and ring with corresponding orbital pattern. B) 
An example of the spin-singlet dimer covering minimizing the 
energy at zero Hund's coupling. Thick (thin) lines denote AF 
(FM) intra- (inter-)dimer bonds, respectively. Dashed lines 
stand for the noninteracting bonds. 



energies E2 = —2 and E^ ~ —3. The dimer states have 
a lower energy than the closed chains if Em > 
where Em is the ground-state energy of the Af-site spin- 
one AF Heisenberg ring. The last inequality is satisfied 
for any AF Heisenberg ring with AI > 4 (see for exam- 
ple Ref. [3]). Note that on a honeycomb lattice only 
the closed chains with M > 6 can be formed. There- 
fore, at zero Hund's coupling 77 = 0, in the ground state 
manifold the pattern of AF bonds corresponds to a hard- 
core dimer covering oi the lattice and on each inter-dimer 
bond only one site has active orbital singly occupied. The 
Haf term is inactive on such inter-dimer bonds, and 
spins of different dimers are decoupled. Spins are cou- 
pled into quantum spin-singlet state on dimer bonds and 
are thus gapped. Such a spin-singlet dimer states form 
the exact ground state manifold of spin-orbital Hamilto- 
nian Eq. 12]) in the limit of zero Hund's coupling. 

Ground state manifold — The ground state manifold is 
extensively degenerate: there are infinitely many ways of 
covering a honeycomb lattice with hard-core dimers, and 
each dimer covering has its own Ising-type degeneracy 
connected with the orientation of an "inactive" orbital. 
Thus, in contrast to the case of one d-electron (S=l/2) 
where the state is determined by the covering of the lat- 
tice with spin-singlet dimers f, here we have an extra 
degeneracy. At each site the second orbital points along 
one of two remaining bonds (the arrows in Fig. IB), with 
the rule that at each of such "empty" bonds there should 
be only one orbital (one arrow). For any dimer covering 
the inter-dimer ( "empty" ) bonds form non- intersecting 
linear chains. Along an inter-dimer path all arrows must 
be directing along the same direction [see as an example 
the zig-zag like linear chain in Fig.[Tj3]. However, on each 
inter-dimer path we can inverse all the arrows and still re- 
main in the ground state manifold. Thus each inter-dimer 
path has a two-fold degeneracy of Ising-type. The degen- 
eracy of a state with given dimer covering is thus equal 



to 2^^^, where Nc is a number of decoupled inter-dimer 



chains. One finds that Nc ^ vN for aU dimer coverings 
except one case when all inter-dimer chains are hexago- 
nal loops. In this case the honeycomb lattice is divided 
into non-overlapping hexagons with no dimers, and the 
remaining hexagons are occupied by three dimers. This 
gives Nc = N/ 6 leading to a finite contribution to a bulk 
entropy from Ising-type degeneracy. The above degener- 
acy is, however, easily lifted by any interaction which in- 
duces a non-zero coupling of non-active orbitals on dimer 
bonds and thus correlates the Ising-like variables (direc- 
tion of arrows) of neighboring inter-dimer chains. For 
example, the coupling of Jahn- Teller distortions on NN 
edge sharing oxygen octahedra favors the ferro-type or- 
bital order of singly occupied non-active orbitals on dimer 
bonds (see Fig. [1)3) [l^. The 7r-type contribution to the 
hopping integral instead stabilizes an antiferro-type or- 
bital order by increasing the AF coupling on a dimer 
bond [20*1. In real materials the two mechanisms will 
have different energy scales and the dominant one will 
dictate the resulting orbital pattern. In both cases one 
completely lifts the degeneracy of a given dimer state. 

The degeneracy of second type originates from the ori- 
entational degeneracy of dimers and is exactly equal to 
the number of hard-core dimer coverings of the honey- 
comb lattice. The latter is well known to be an extensive 
quantity. The classical dimers on a honeycomb lattice ex- 
hibit a power-law decay of dimer-dimer correlation func- 
tion and are thus in a critical state j2l|. Therefore, one 
expects that any perturbation which may favor one or 
another type of dimer orientation may lift this extensive 
degeneracy, resulting in a long-range ordered pattern of 
dimers. We have considered a possibility of perturbing 
the exact dimer ground state manifold by small but finite 
Hund's coupling < ry <C 1. The Hund's coupling pro- 
duces a ferromagnetic (FM) term i/pM active on inter- 
dimer bonds (see Fig. [Tj3). At small ry, the dimer state is 
stable against the weak FM inter-dimer interaction. In 
this case the magnetic contribution along the FM bond 
is zero {{Si ■ Sj) = for i and j belonging to different 
dimers). However, the weak inter-dimer coupling may, in 
principle, lift the orientational degeneracy through order 
out of disorder by triplet fluctuations [o'] . We find that on 
a honeycomb lattice these quantum fluctuations do not 
fully lift the dimer degeneracy. They only disfavor the 
configurations with three dimers on a hexagon because 
such a hexagonal loop of alternating AF and FM bonds 
is frustrated in the classical limit. In the next section 
we introduce a magnetoelastic coupling and show that it 
fully lifts the orientational degeneracy of dimers. 

Degeneracy breaking by magnetoelastic coupling - The 
magnetoelastic mechanism of lifting the extensive degen- 
eracy of frustrated systems has been successful in explain- 
ing the experimentally observed structures (I,!!!, 23. 2^. 
The physics behind this mechanism is that the correlated 
nature of structural distortions that appear due to the 





FIG. 2: (Color online) Lifting dimer degeneracy by magne- 
toelastic coupling. The small light (large dark) circles denote 
Ru (Li) ions. The arrows show the displacements of corre- 
sponding ions. A) The sketch of Li displacements induced by 
dimerized Ru-Ru bond. B) and C) Two examples of destruc- 
tive interference of Li displacements induced by neighboring 
dimers. D) The ground state dimer pattern selected by mag- 
netoelastic coupling. This pattern exactly corresponds to the 
one found in Li2Ru03 in Ref. 0). 



modulation of magnetic energies on bonds may select a 
particular pattern of distortions and thus lift the exten- 
sive degeneracy of the ground state manifold. 

In the ground state manifold the spin degrees of free- 
dom are described by the product of spin-singlet dimer 
states. A shortening of a bond where the singlet is lo- 
cated enlarges the magnetic energy gain, because of the 
increase in the exchange coupling J on that bond. The 
magnetic energy gain, being linear in distortion, out- 
weights the increase in elastic energy and would always 
lead to such a contraction of singlet bonds for any dimer 
covering. But due to elastic coupling of these distortions 
different distorted patterns will lead to different elastic 
energy and hence to the lifting of dimer degeneracy. We, 
first, formalize this picture within the simplest model and 
later discuss its extension. The schematic structure of 
Li2Ru03 is shown in Fig. 2D. In the minimal model the 
exchange coupling is assumed to depend solely on a dis- 
tance between Ru ions, and the lattice degrees of freedom 
are described within Einstein phonon model for Ru ions: 



ME 



2 ^ ' 



(3) 



where 7 = " 7 ^'^dr^ ^ magnetoelastic coupling con- 
stant, e is magnetic energy gain on a dimer bond, Ui is 
a displacement of Ru ions at site i, and k is Einstein 
phonon constant. The dimer positions are described by 
lij: the latter is an unit vector along a bond ij occu- 
pied by dimer and is zero otherwise. We minimize the 
energy functional £me with respect to set of Ui and elim- 
inate them in favor of lij . We find £me — ~^ Si '^i > ^'^'^ 
Ui = TP X] 7 • The set of Ui define the distortion pattern 
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in which Ru ions on dimer bonds are moved towards each 
other. However, as the dimer bonds do not share a com- 
mon site, ^ME is a constant. It is independent of dimer 
variables . In order to introduce the couphng between 
the dimers we extend the model by including a finite force 
induced on Li ions by distorted pattern of Ru ions. Con- 
sider an isolated complex of a Ru-Ru bond together with 
neighboring Li ions, shown in Fig. [2]A_. When this bond 
is occupied by spin-singlet dimer it gets contracted, and 
the displacements of Ru ions induce the corresponding 
distortion of Li complex as seen in Fig. O'V. On a lattice 
the distorted dimer bonds will now be coupled through 
the force induced on a common Li ion, an effect similar to 
a cooperative Jahn- Teller physics. The dimer configura- 
tions for which the induced forces on Li sites interfere in a 
non-destructive manner result in a larger distortion and 
hence more gain in energy. The two orientations of neigh- 
boring dimers shown in Fig. [33 and C are energetically 
unfavorable, as in both cases the induced forces exactly 
cancel each other. Thus the ground state dimer pattern 
should satisfy the constraint of noB and no-C type con- 
figurations. Note that a dimer covering satisfying the 
former constraint automatically satisfies the latter. It is 
possible to check that the only possibility to fulfill such a 
no--B constraint is realized for the dimer pattern shown 
in Fig. [2jD. This dimerization pattern exactly reproduces 
the one observed in the insulating phase of Li2Ru03 Q. 

We note that above described mechanism of the selec- 
tion of spin-singlet dimer pattern equally applies to the 
case when singlets are formed by spin-one-half degrees 
of freedom. The spin-singlet dimer nature of the ground 
state manifold for systems on a honeycomb lattice 
has been proven in Ref. We therefore predict the 
same dimerized superstructure also for systems, such 
as V*"*" and Ti'^+ based compounds with such structure 
(if orbital degeneracy will not be lifted by some external 
mechanism like trigonal distortion, often present in such 
structure). 

Summary.- To summarize, we have studied the ground 
state of spin-one honeycomb antiferromagnet with par- 
tially filled levels. We have demonstrated that the 
orbital degeneracy induces spontaneous dimerization of 
spins and drives them into extensively degenerate man- 
ifold of spin-singlet dimer states. The orientational de- 
generacy of dimers is then lifted through the magnetoe- 
lastic interaction that stabilizes a peculiar valence bond 
crystal state. Our theory provides an explanation for 
the observed nonmagnetic dimerized superstructure in 
Li2Ru03 compound. 
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